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Abstract. 

We study the arrival directions of 69 ultra-high energy cosmic rays (UHECRs) observed 
at the Pierre Auger Observatory (PAO) with energies exceeding 55 EeV. We investigate 
whether the UHECRs exhibit the anisotropy signal expected if the primary particles are 
protons that originate in galaxies in the local universe, or in sources correlated with these 
galaxies. We cross-correlate the UHECR angular positions with the positions of IRAS-PSCz 
and 2MASS-6dF galaxies (with median depth of 120 and 225 Mpc respectively) weighted for 
GZK and other particle propagation effects. This is the first time that the 6dF survey is 
used for a study of this type and the first time that the PSCz survey is used with the full 69 
publicly released PAO events. We find the probability that the observed UHECR events are 
consistent with an isotropic hypothesis to be less than 5% for most of the parameter space 
considered. The observed correlation is consistent with the mean of the model distribution 
of UHECRs associated with PSCz and 6dF galaxies. The agreement in results between the 
two catalogues, which probe similar large scale structure, is very good. Finally, we explore 
how the random magnetic deflections of UHECRs during propagation affect the expected 
anisotropy signal. 
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1 Introduction 

The origin of cosmic rays (CRs) with energies exceeding 10 EeV (1 EeV = 10 18 eV) remains 
unknown despite decades of extensive research. Such particles are referred to as ultra-high 
energy cosmic rays (UHECRs) and their sources are probably extra-galactic [1], [2]. Their 
gyro-radius is larger than could be contained in the magnetic field of our galaxy [3]. 

UHECRs with energy above ~ 50 EeV are above the threshold for pion photo-production 
upon collision with cosmic microwave background (CMB) photons, a process which was 
predicted soon after the discovery of CMB radiation, known as the GZK cutoff [4], [5]. The 
most recent published cosmic ray spectra measured by the High Resolution Fly's Eye (HiRes), 
PAO and Telescope Array (TA) experiments confirm this cutoff in the CR spectrum [6] , [7] , 
[8]. As a result of the GZK process UHECRs that arrive on earth with energy exceeding 
50 EeV must originate in sources within a few hundred mega-parsecs (Mpc). 

Only a handful of known astrophysical populations are likely to have the required power 
to accelerate CRs to such high energies [3] , [9] . The most likely candidates are active galactic 
nuclei (AGN) and gamma ray bursts (GRBs) (see however [10], [11], [12] for a recent result 
and alternative interpretations on GRBs as UHECR sources). 

If extra-galactic magnetic fields (EGMFs) are not too large, UHECRs must point back 
to their sources within a few degrees if they are protons, due to their high rigidity. Hence if 
they originate in some astrophysical population, the arrival direction distribution of UHECRs 
should be correlated with that population as well as with the distribution of large scale 
structure (LSS) in the local universe, since matter in the Universe is clustered [13]. An 
isotropic arrival direction distribution of UHECRs could indicate either a heavier composition 
or new physics such as heavy relic particles (see [1] for a review of experimental signatures). 
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Besides the origin, another question which remains open is the composition of UHECRs. 
The PAO, the largest UHECR observatory to date, has not reported on the composition of 
UHECRs with energy exceeding 20 EeV, presumably due to the small number of observed 
events. Their most recent results between 1 — 20 EeV suggest that there is a smooth transition 
from light, proton-like to heavier, Fe-like nuclei [14]. It should be noted that HiRes and TA 
results up to an energy of 50 EeV agree with a proton-like composition [15], [16]. There 
are a number of reasons why one wouldn't expect heavy nuclei at energies above the GZK 
threshold (see for example [17]) not least that the leading extra-galactic candidate sources 
are expected to accelerate primarily protons. 

In this work we revisit the question of the origin of UHECRs, following the release of the 
arrival directions of 69 UHECRs with energy above 55 EeV recorded up until December 2009 
by the PAO [18]. We model the local UHECR source distribution using galaxy catalogues of 
the nearby universe. The 6 degree Field Galaxy survey (6dF) is being used in this paper for 
the first time to derive the expected UHECR source distribution. The IRAS Point Source 
Catalogue of redshifts (PSCz), although previously used for cross correlations with UHECRs 
(see references below), is being used here for the first time to analyse the updated dataset 
of 69 observed UHECRs. We refrain from imposing any artificial distance cuts on potential 
UHECR sources and instead introduce a weighting system which ensures that the expected 
contribution of distant sources to the UHECR flux on earth drops off naturally with distance. 
We also investigate the effect of random magnetic deflections on any correlation between the 
source population and the observed UHECR arrival direction distribution. 

A number of authors have looked for the sources of UHECRs in the past, by cross- 
correlating the observed arrival directions with the positions of nearby galaxy catalogues, 
in particular the PSCz [19], [20], [21], [22], [23], 2MRS [24], [18], [25] and others with the 
positions of specific classes of objects, AGN [26], [27], [28], [29], [30], BL LAcs [31] and 
LIRGs [20] among others, reporting different degrees of correlation depending on the UHECR 
sample, statistical approach and source population used. The authors of [19] analysed the 
arrival directions of 27 UHECRs with energy above 55 EeV detected until August 2007 and 
found that the UHECRs exhibited a stronger correlation with propagation weighted PSCz 
galaxies than 99.8% of isotropic realisations. The authors of [24] analysed the distribution 
of the 69 observed UHECRs and found that 21 of the 55 that survive their cuts correlate 
with nearby AGN. The probability of finding such a correlation assuming isotropy is 0.003. 
Further they cross-correlated the arrival directions of the UHECRs with the positions of 
2MRS galaxies and Swift-BAT X-ray sources [32]. For the values of the free parameters in 
their models that maximise the likelihood they found that the fraction of isotropic realisations 
that yield a higher likelihood than the observed UHECRs are 0.004 and 2 x 10~ 4 for the 2MRS 
and Swift-BAT respectively. The results they obtained are a posteriori and do not constitute 
a confidence level on anisotropy. Because of the different approaches followed their results 
are not directly comparable to ours. 

This work is organised as follows. In Section 2 we present the formalism used to analyse 
the arrival direction distribution of UHECRs, in Section 3 we present our main results applied 
to the measured PAO data and in Section 4 we discuss the implications of our results and 
conclude. 
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Figure 1: The function P p (vl, Ef, Ef), which represents the probability that a UHECR 
emitted at a source at luminosity distance ri, with initial energy Ei, will be observed with 
final energy above Ef, here shown as a function of propagation distance for Ei = 200 EeV 
and E f > 40, 60, 80, 100 EeV. 

2 Methodology 

We study the correlation between the arrival directions of the 69 PAO observed events (§2.1) 
and the distribution of matter in the local LSS which we model using the 6dF (§2.2) and 
the PSCz (§2.3) catalogues. In §2.4 we present our models of UHECR source distribution. 
Throughout this analysis we take into account the non-uniform PAO exposure (§2.1). In §2.5 
we explain our treatment of the expected UHECR random magnetic deflections. In §2.6 we 
present our statistical approach. 

2.1 The Pierre Auger observatory 

The PAO, in Malargue Argentina is a cosmic ray observatory dedicated to the detection of 
cosmic rays of energy greater than of 1 EeV. Data taking started in January 2004 and since 
then the arrival direction and reconstructed energy of 69 UHECRs exceeding 55 EeV have 
been published. At the PAO, which has a detection area of 3000 km 2 , UHECR particle show- 
ers are detected by 1600 ground based Cherenkov detectors surrounded by 24 fluorescence 
telescopes, which measure the amount of energy dissipated in the atmosphere in the form 
of ultra-violet radiation [33]. The combination of these two techniques provides the most 
accurate reconstruction of UHECR shower geometry to date. Full time operation ensures 
an exposure which is uniform in right ascension and fully efficient for zenith distance up to 
9 m = 60° . The PAO exposure is a function of declination and is given by [34] : 



lo(5) oc cos(ao) cos(<5) sin(a m ) + a m sin(ao) sin(<5) 



(2.1) 
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where 5 is the declination and ao is the latitude of the PAO which is —35.2°. Here, a m is 
given by: 

'0 if£>l 
a m = < vr if £ < -1 

i cos~ 1 (^) otherwise 



and 



cos 



sin(ao) sin(<5) 



cos(ao) cos(<5) 



The PAO has an integrated exposure of ~ 20000 km 2 year and ~ 20 UHECRs above the 
GZK threshold are detected every year since its completion in November 2008 [35] . 



2.2 The 6 degree field galaxy survey 

The 6dF is a redshift and peculiar velocity survey of 2MASS selected galaxies which was 
carried out using the Six-Degree Field instrument on the Schmidt Telescope of the Anglo- 
Australian Observatory [36] . The survey, which covers the entire Southern Sky (excluding the 
Galactic plane, \b\ > 10°), resulted in a catalogue of 125,071 galaxies and 110,256 associated 
redshifts. The 6dF has a median redshift z = 0.053 which corresponds to a comoving distance 
of 225 Mpc for a flat universe with Qm = 0.25, = 0.75 and H Q = 70 km s _1 Mpc -1 , which 
we assume throughout this work. The 6dF field of view covers 80% of the PAO field of view 
by area. Taking into account the total declination dependent PAO acceptance (equation 2.1) 
which is smaller for positive declinations, it covers 86.2% of the instantaneous PAO exposure. 
The 6dF is near-infrared selected, which means it favours older, bulge-dominated galaxies 
and therefore it is minimally affected by dust extinction. In this work the final 6dF data 
release [36] is being used for the first time for a study of the origin of high energy particles. 



2.3 The IRAS PSCz 

The PSCz is a redshift catalogue of infrared selected IRAS galaxies which covers 84% of the 
sky [37]. It contains 14,677 galaxies with associated redshifts, and median redshift (z = 0.028) 
which corresponds to a comoving distance of 119 Mpc in the cosmological model we are 
assuming. Since the galaxies in the PSCz are infrared selected there is a preference for 
young, star-forming galaxies in the catalogue. 

We have chosen to use two complementary galaxy surveys to derive the expected 
UHECR source distribution. The PSCz is a shallow nearly full-sky galaxy survey and its 
use facilitates comparison with results of previous studies. The 6dF on the other hand is a 
much larger survey (~ 20 times more galaxies than the PSCz in the southern hemisphere). 
The different median depths of the 2 surveys mean that they highlight different structures of 
the nearby universe (for example the Shapley Concentration (centred at (I ~ —50°, b ~ 30°) 
at a distance ~ 200 Mpc is prominent in the 6dF) although there is a significant over- 
lap. The 2 surveys consist of different galaxy populations with slightly different clustering 
properties. This means that they can be used to distinguish between different astrophysical 
populations as UHECR sources, although this is a subtle difference that will require a large 
UHECR dataset and a better understanding of UHECR magnetic deflections before it can 
be pursued. 
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2.4 Model of UHECR source distribution 

We use the galaxy catalogues introduced above to model a UHECR source population which 
is steady and follows the distribution of luminous matter in the local universe. One of the 
parameters that our analysis is sensitive to is the density of UHECR sources. The uncertainty 
on the true position of the sources of observed UHECRs due to random magnetic deflections, 
means that groups of two or more UHECRs with arrival directions separated by less than 
a few degrees may have been produced by the same source and are referred to as repeaters. 
The number of repeaters in a dataset can constrain the density of sources of UHECRs [13] , 
[38], [39]. We have searched for repeaters separated by less than 3°, intended to reflect the 
uncertainty associated with instrumental resolution (~ 1°) and possible magnetic deflections 
of a few degrees (see §2.5), among the 69 observed PAO events. We found 9 doublets and 2 
triplets, which suggests ~ 600 UHECR sources within the earth's GZK horizon and hence a 
source density ~ 4 • 10 -4 Mpc~ 3 . 

Here we consider a model in which the local number density of UHECR sources is 
comparable to that of galaxies in the local Universe (no = 10 -2 Mpc~ 3 ) and UHECR sources 
are faint i.e. each source produces 1 or no events and the probability of a single source 
producing multiple events is low. If only one subtype of galaxy is responsible for the observed 
UHECR flux or the observed flux is dominated by a few "bright" sources (such as in models in 
which the nearby radio galaxy Centaurus A is responsible for much of the observed UHECR 
flux [40], [41], [42]) then the density of sources may be lower than considered here. More data 
will allow us to better constrain the UHECR source density. In our model we also assume 
that all UHECR sources are intrinsically identical. 

We assign each galaxy in our model a weight, proportional to the expected UHECR flux 
from that source, which depends on the flux suppression with distance and the energy with 
which protons are emitted (this initial energy determines the energy losses during propagation 
due to interaction with background photons). The flux weight for a galaxy at luminosity 
distance from earth that emits a cosmic ray proton with initial energy E% which reaches 
the earth with energy equal to or greater than Ef, is given by: 



is the derivative of the function Pp{ri, Ef, Ef) which gives the probability of a proton arriving 
on earth with energy Ef if it was emitted with energy Ei by a source at distance tl, introduced 
in [43]. For Ef we conservatively adopt the lowest measured energy (55 EeV) present in the 
PAO dataset. We set E L max which is the maximum energy achievable through astrophysical 
processes equal to 10 21 eV. The value of E^max is not fully constrained by observations but 
considering higher values of -E^max would have a negligible effect on our results here, since 
UHECRs above this energy promptly interact with background photons. For the intrinsic 
spectrum of the UHECR sources I(Ei,t) we consider a power law spectrum: 




(2.2) 



where 



p p (r L ,Ei;Ef) 



dP p (r L ,E i; E f ) 
dE f 



(2.3) 



I{Ei) = hEf@{E- K 



i,max 



Ei). 



(2.4) 
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The step function restricts the UHECRs to have energy below -Ei jmax - Here Iq normalises the 
injection spectrum such that: 

E- 

dEiI E?e(E iimax - Ei) = 1. (2.5) 

We consider a power law with index a = — 2 which is what is expected in diffusive shock 
acceleration [44] , [45] . We have investigated the sensitivity of our model to a and found that 
varying the value of a in the range —2.5 < a < —1.5 has a negligible effect on our results. 
For the simulation of the proton propagation and energy losses we use the |»(r_., Ef, Ef) 
function which was computed and made publicly available by the authors of [39] and [46] 
at http://www.desy.de/~uhecr/P_proton. In figure 1 we plot the function P p (tl, Ef, Ef), 
which we obtain by numerically integrating over the p(rL,Ef, Ef) function, for UHECRs with 
Ei = 200 EeV and E f > 40, 60, 80, 100 EeV. 

To account for the effect of the flux limit of the survey on the observed number density 
of galaxies as a function of distance we need the survey's selection function, defined as 
the expected number density of galaxies in the survey as a function of distance, in the 
absence of clustering. One can model the selection function using a fit to the survey's redshift 
distribution, often parametrised as: 



dN(z) = Az p 



exp 



> 7 
Z x 



dz, (2.6) 



where A gives the normalisation, z p the peak of the distribution and f3, 7 control the slope. 
All 4 parameters are specific to the survey (see [36] and [37] for 6dF and PSCz redshift 
distributions respectively). The overall selection function ip(r c ) is the redshift distribution 
divided by the volume element [47]: 

. . 1 [dN\ ( dz\ 



Q. s rl \dz J \dr, 

where £l s is the solid angle of the survey and r c the comoving distance. We weight the 
observed galaxy distribution by the inverse of ip(r c ) which we normalise to be 1 at r m j n = 
5 Mpc comoving, so that the effective contribution of each survey galaxy to the model source 
distribution is: 

_ u(r L ) flux 
ip[r c ) 

Galaxy surveys measure redshifts not distances, and peculiar velocities along the line- 
of-sight affect our estimates of actual distances. The redshift of an object is defined by: 

cz = H r + (v (r) - v (0)) • f , (2.9) 

where v (r) is the object's peculiar velocity and v (0) is the observer's peculiar velocity. 
Working in a reference frame where Av = v (r) — v (0) is small allows more accurate distance 
estimates. In the local universe, out to cz ~ 3000 kms -1 , where galaxies share the motions 
of the Local Group (LG) it is best to convert to the LG rest frame. Further away galaxy 
peculiar velocities are independent of the LG velocity and Av is smaller in the CMB rest 
frame. In this work we are interested in the distribution of matter nearby, hence it is most 
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useful to work with LG frame redshifts. However we will also present our main results from 
redshifts in both frames. 

To convert between czuei, that is redshift measured in the heliocentric frame, and czlg 
we use the relation [48]: 

cz LG = cznei - 79cos(/) cos(6) + 296sin(/)cos(6) - 36sin(&), (2.10) 

where I and b are the Galactic longitude and latitude respectively. To convert from the LG 
frame to the CMB frame we use: 

czcmb = cz LG + v LG [sin(fc) sin(&L G ) + cos(6) cos(b LG ) cos(\1 L g ~ , (2.11) 

where vlq is the amplitude of the LG velocity with respect to the CMB and {IlGi^lg) is the 
direction of its motion. Using the CMB dipole value obtained in [49] with 7 year data from 
WMAP we obtain v LG = 627 km s" 1 towards (l LG = 276°, b LG = 30°). 

2.5 Magnetic fields 

A relativistic particle of energy E and charge Z in a magnetic field has a Larmor radius 
given by ri ar = z f B , where B± is the field strength in the direction perpendicular to the 
momentum of the particle. A relativistic proton in a field which is constant over a distance 
d will suffer a deflection 

{^Y\ik)(whY (2 - 12) 

The strength and distribution of EGMFs are poorly known (see reviews of existing data 
in [50], [51]). Dense large scale structures such as galaxy clusters and galaxy filaments are 
likely to support relatively large magnetic fields whereas in between structures EGMFs are 
probably negligible. The deflection for a proton propagating a distance d through a magnetic 
field with r.m.s. strength B and coherence length A is given by [52]: 

The most realistic models of UHECR propagation through EGMFs model deflections as a 
collection of scattering events in random directions. Large scale simulations of EGMFs [53], 
as well as results from a recent semi-analytic analysis [54] , conclude that deflections of proton 
UHECRs of energy E > 40 EeV do not exceed 3° over 99% of the sky for d ~ 100 Mpc. 

Magnetic deflection which is proportional to the charge is much larger for heavier 
UHECR nuclei and can completely wash out the directional correlation of UHECRs with 
their sources [55]. The composition of the UHECR sample we are considering here is at 
present uncertain but the anisotropies we would expect at lower energies if this dataset was 
composed of nuclei rather than nucleons have not been observed (see [56], [57]). Through- 
out most of this analysis we consider proton UHECRs. In order to bracket possible proton 
UHECR magnetic deflections we perform our analysis by averaging over angular bins in the 
range 3.9° — 7.3°. We exclude the Galactic plane fom our analysis, \b\ > 12°, since UHECRs 
travelling through the Galactic plane may suffer strong deflections due to the magnetic field 
in the Galactic disc. In §3.2 we relax the assumption of nucleonic composition and inves- 
tigate the effect that larger random deflections have on any correlation between the source 
population and the observed arrival direction distribution. 
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2.6 Correlation statistic 



To detect any existing anisotropy signal we consider the counts-in-cells statistic X, proposed 
in [19], which characterises the correlation between the predicted and observed UHECR 
arrival direction distribution. In [19] X was shown to be more sensitive to the expected 
anisotropy signal than the angular power spectrum (e.g. [34], [58]) and the two point correla- 
tion function [59], which are other statistical measures commonly used in clustering analyses. 

Let WAuger.i be the exposure of PAO in bin i (see equation 2.1) and CL) survcyi i the weight 
imposed by the survey mask (given by [36], [37]): 



^6dF,i 

WpSCz,i 

We define a combined weight: 



1 for \b\ > 10° and 5 < 0° 
otherwise 

fl for |6| > 10° 
I otherwise 



^exposure, i — ^ Auger, i ' ^survey, i- (2T4) 

We require that each of the quantities that go into our statistic are (a) weighted exactly once 
by the combined weight and (b) normalised so that each has a sum equal to the weighted 
sum of cosmic rays in the mask defined region. We restrict the analysis to the mask defined 
region by subjecting the number of cosmic rays in each bin i in the PAO dataset AA U ger,i to 
the survey mask: 

NcR,i = W surveyj i • iVAuger.i- (2-15) 

Note that iVAuger,ii the observed number of cosmic rays, is naturally subject to WAuger.i so we 
don't apply this weight again. 

The number of survey galaxies visible by the PAO, is limited by the observatory's 

declination dependent exposure (equation 2.1). To account for this effect we subject each 
galaxy in bin i to the combined weight: 

NgaLi — W ex p 0sure i • \ — — : r— , (2.16) 

ip(r c )j 

where the sum is over the weighted contribution to the flux of each survey galaxy j in bin 
i. In the model where UHECR sources are correlated with galaxies in the nearby LSS, the 
expected number of cosmic rays iV~M,i is simply N gaj i,i normalised to the number of observed 
UHECRs according to condition (b): 

AT M ,i = • AW (2.17) 

In an isotropic model the expected number of cosmic rays is given by: 

= ^5H£ 2A . ( ^ -E^AugerJ ] , (2.18) 

where n is the number of angular bins in the sky. The integral J 2ir s'm(5)w(5) d5 gives the 
total exposure of the PAO and hence j^ 27 rsin(S)w(8) d5 " unmas ks" A^Auger.i, i.e. the quantity 
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in the brackets gives the number of cosmic rays that would have been observed if the PAO 
had had a uniform, full sky acceptance. 
The statistic X is defined as: 



We generate mock realisations of isotropic UHECR samples and correlate them with the 
galaxy maps to model the behaviour of X in an isotropic model. Similarly we create mock 
realisations of UHECRs that trace the local LSS by sampling from a Poisson distribution 
with a mean equal to the number of UHECRs expected to be detected per unit area at the 
PAO above energy Ef using the model of UHECR source distribution presented in §2.4. 
Finally we compute the statistic for the observed PAO UHECRs, and determine which of the 
2 models of source distribution is preferred by the observed PAO UHECR dataset. 



The methodology presented in Section 2 was applied to the 69 UHECRs observed by the PAO 
with energy E > 55 EeV until the end of 2009 [18]. We perform our analysis by dividing 
the sky in equal area bins (counts-in-cells) . The disadvantages of the counts-in-cells scheme 
are the effect of boundaries and the ability to arbitrarily choose the bin size. We deal with 
both these limitations by treating the size of the equal area bins as a free parameter in a 
range that covers expected random proton UHECR magnetic deflections (§3.1). In §3.2 we 
relax the assumption of proton UHECRs and investigate the sensitivity of our results to the 
magnitude of the magnetic deflection. We note that as explained in Section 2 we always 
mask out the galactic plane, as well the regions of the sky not covered by the galaxy survey, 
so the number of observed UHECRs that remain in the mask defined region is never 69 and 
depends on the survey used and the bin size considered. For simplicity we still write 69 
events to specify the dataset that is being used. 

3.1 Cross-correlation of UHECRs and nearby LSS 

In figure 2 we show the distribution of sources of 55 EeV UHECR which follow the distri- 
bution of matter in the PSCz/6dF catalogues. To aid with illustration, we plot the model 
source distribution smeared with a 2 dimensional Gaussian filter with standard deviation 
a = 3°, 7.2°, whereas throughout the rest of this work we have treated galaxies as discrete 
points. We plot the source distribution with the 69 observed PAO UHECRs superimposed in 
figure 3. In the maps of the PSCz one can see a very large contribution from the Virgo cluster 

(Z 80°, b ~ 75°), Hydra-Centaurus (-60° < I < 0°,0° < b < 45°) and the Perseus-Pisces 

supercluster (I ~ 140°, b ~ —25°). In the maps of the 6dF there is a large excess as a result of 
flux from Hydra-Centaurus and the Shapley Concentration (centred at (I ~ —50°, b ~ 30°) at 
a distance ~ 200 Mpc). Note that we haven't masked with the PAO exposure in these plots. 
After doing so the relative weights of structures change so that for example Virgo no longer 
dominates in the PSCz derived model. Sky maps have been created using the HEALPix 
package [60]. 

Because the observed galaxy distribution is in the form of flux limited as opposed to 
volume limited samples, which would be more suitable for our analysis, we have weighted 
observed galaxies in the PSCz and 6dF with the inverse of the selection function ip{r c ) as 




(2.19) 



3 Results 
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0.064 3.1 0.0 3.4 



Figure 2: Maps of the predicted, local UHECR source density, in Galactic coordinates 
centred at the Galactic centre with the longitude / increasing anti-clockwise, derived from 
the PSCz (left) and 6dF (right) for UHECRs with final energy 55 EeV. The predicted source 
distribution has been smoothed with a Gaussian filter, with a = 3° (top row) and a = 7.2° 
(bottom row) for presentation purposes. The intensity at each point is normalised to its all 
sky average (half-sky for the 6dF). The PAO is sensitive to the part of the sky below the 
thick black line. 




Figure 3: The predicted UHECR source density as in figure 2 above, but with the 69 PAO 
observed events superimposed (in grey). The PAO exposure has not been taken into account 
in this plot, hence visual inspection of correlations can be misleading. A statistical analysis 
as in §2.6 is required. 



outlined in §2.4. There is a degree of uncertainty associated with this process, i.e. an error 
associated with the selection function which we have not taken into account, however the 
resulting flux contribution is likely to be close to its true value since matter in the universe 
is clustered. 

In figures 4 and 5 we show the distribution of values that the statistic X takes in 
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Figure 4: The distribution of values of the correlation statistic X (equation (2.19)) in 
10,000 mock realisations of a set of 69 UHECRs drawn from an isotropic distribution (dashed 
histograms) and from the predicted UHECR source distribution generated using the PSCz 
(solid histograms). The black vertical line indicates the value of X obtained for the 69 
observed PAO events (^Auger)- Each subplot corresponds to a different cell size used for the 
counts-in-cells analysis in the range 7.3° (top left) to 3.9° (bottom right). 
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Figure 5: The same distributions as in figure 4 but using the 6dF catalogue. 



10,000 mock realisations of 69 isotropically distributed UHECRs, as observed by the PAO 
(dashed histograms), and the distribution of values of X in 10,000 mock realisations of 55 
EeV UHECR protons with sources drawn from the predicted UHECR source distribution in 
the case where UHECR sources follow the distribution of matter in the PSCz and 6dF cata- 
logues respectively (solid histograms). The black vertical line in each subplot shows the value 
of X for the 69 observed PAO UHECRs assuming Ef = 55 EeV. Each subplot corresponds 
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Figure 6: Left: Pi so - dark grey (Ppscz - light grey), the percentage of realisations of an 
isotropic (correlated with LSS) UHECR source distribution in which the value of X was 
higher than the value of X obtained for the 69 observed PAO UHECRs, as a function of 
cell-size. Right: Same as on the left but using the 6dF survey. Results shown in solid lines 
are obtained using redshifts in the LG frame and in dashed lines with redshifts in the CMB 
frame. The horizontal dotted line shows the 95% confidence level. 
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Figure 7: Ratio of the frequency with which the value of X obtained with the observed 
PAO UHECRs (AAuger) is obtained in each of the two models of UHECR source distribution 
we are considering (following LSS -LSS, isotropic -ISO) p^lojpAoj as a function of cell-size, 
for a model source distribution that follows the PSCz (left) and 6dF (right). 



to a different cell-size in the range 7.3° — 3.9°. We see that throughout the parameter space 
considered the observed PAO UHECRs are consistent with a source distribution correlated 
with nearby LSS whereas they have a value of X higher than in most mock isotropic UHECR 
realisations. We also see good agreement between the results obtained with the two galaxy 
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Figure 9: Distributions of the statistic X for realisations of UHECR sets with isotropic 
arrival directions (dashed histograms) and arrival directions correlated with the LSS (solid 
grey histograms). The value of X for the 69 observed PAO events is given by the black 
vertical line, for the 6dF (left), PSCz (centre), and PSCz only in the southern hemisphere 
(right). The plots shown are for cell size 5.9° x 5.9°. 



surveys. Although the 6dF and PSCz samples are selected differently (near-infrared and 
far-infrared respectively) and the 6dF is deeper than the PSCz, they probe approximately 
the same large scale structure. 

We quantify the significance of this result in figure 6, where we show the probability 
that a value of X greater than that obtained for the observed PAO UHECRs (^Auger) will 
be obtained in each of the models of UHECR source distribution (Pi so in the isotropic model 
and -PpsCz/-P6dF hi the model where the UHECR sources are correlated with the local LSS), 
as a function of cell-size. Throughout the parameter space considered the observed PAO 
UHECRs exhibit a stronger correlation with nearby LSS than 90% of all isotropic mock real- 
isations and are consistent with a source distribution that follows the distribution of matter 
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in the PSCz/6dF surveys. In the same figure we show the results obtained using redshifts 
in the CMB as well as the LG frame. There is a larger spread between the predictions from 
the two reference frames for the PSCz catalogue, which is as expected, because the flux is 
dominated by more nearby structures. 

There are differences between the model source distribution of ref. [19] derived from 
the PSCz catalogue and the PSCz derived model presented here, a main one being that they 
have considered a smoothed galaxy distribution whereas we have considered a distribution of 
point sources. When analysing the original 27 observed UHECRs using the PSCz to model 
the source distribution with comparable parameters, Ef = 55 EeV (57 EeV in [19]) and cell 
size 6.3° x 6.3° (6° x 6° in [19]), we find P iso = 1.03% and P PS Cz = 8.7%, which are similar 
to their P 7 = 0.9% and P UB = 3.8%. 

Another way to quantify the significance of the results of figures 4 and 5 is to con- 
sider the frequencies with which the value of X obtained for the observed PAO UHECRs 
occurs in mock realisations of the two different models of UHECR source distribution. If 
the ratio of the likelihoods of the two model source distributions given the observed out- 
come p(Modcl ^outcome) * s S rea ter than 1 then Model 1 is preferred by the observed data. We 

i , i ■ ,. P(LSS source model is true|JfAu E cr) 1 c±. P(LSSIPAO) , i r ,1 , 

show this ratio P(isotropic source modcl is truc |xL gcr ) » hereafter P | Iso j PAO) ; separately for the two 
galaxy catalogues as a function of cell size in figure 7. The results shown in this plot vary 
considerably with cell size, but throughout the parameter space we have investigated there 
is a preference for a distribution of UHECR sources correlated with LSS, consistent with the 
results of figure 6. 

Despite our treatment of our own peculiar velocity in space, there remain galaxies in 
the surveys we use to model the UHECR source distribution which appear to be much closer 
to the earth than they are in reality and others that appear to have a negative redshift as a 
result of their own peculiar motions. To deal with this we assign galaxies whose distances are 
known from independent measurements their measured distances, instead of converting them 
from redshifts and we remove from the model of UHECR source distribution all galaxies that 
are (or appear to be) within less than 5 Mpc comoving. In order to test the robustness of our 
results to this cut we consider how cuts of 10 and 15 Mpc affect our results. In figure 8 we 
plot the quantities Pi so and -PpsCz/-P6dF (defined earlier on in this section) after removing all 
galaxies within 5, 10 and 15 Mpc form the UHECR source model as a function of cell size. 
The effect of this cut is almost negligible in the 6dF analysis because it is a slightly deeper 
survey and flux from the nearest galaxies does not dominate the expected flux on earth. The 
results obtained with the PSCz are more sensitive to this cut, but our conclusions are robust 
with respect to this cut since in all cases and for both surveys the PAO observed UHECRs 
have a higher value of X than 90% of all isotropic mock realisations and are consistent with 
a distribution of sources which follows LSS. 

In figure 9 we show the results of the correlation analysis performed using the 6dF sur- 
vey next to the same results for the PSCz, as well as the results of the analysis performed 
with the PSCz survey but restricted to the Southern hemisphere. This allows for a direct 
comparison between the results of the two surveys since we confine the analysis to the part 
of the sky that they are both sensitive to. Again we see good agreement between the two 
surveys and the observed PAO UHECRs are consistent with a distribution of sources corre- 
lated with nearby galaxies. 

There is a 23% systematic uncertainty in the energy determination of the primary par- 
ticles detected at PAO which introduces an error to our results. We have investigated how 
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Figure 10: The distribution of values of X in 10,000 realisations of 69 isotropic UHECR 
events (dashed histogram) and 69 UHECRs whose sources are in local galaxies (solid his- 
tograms) with simulated magnetic deflections (PSCz- Ze/t subplot and 6dF-right subplot). In 
each solid histogram the UHECRs arrive on Earth with a mean position centred at the 
source's true position and random angular deflections with amplitude a given in the legend. 

this uncertainty affects our results and found that they are robust to this uncertainty, i.e. 
changing Ef in our models by 23% does not change our conclusions. 

3.2 Magnetic deflection 

We investigate the dependence of our results on the amplitude of the magnetic deflections 
suffered by UHECRs during their propagation. If UHECRs are Fe nuclei as some experi- 
mental data suggest [18], average deflections are likely to be of the order of 20°. The results 
obtained in this section are relevant to protons and Fe nuclei which are attenuated over sim- 
ilar distances at energy above ~ 40 EeV but not to intermediate mass nuclei which have a 
much shorter mean free path at this energy range (see for example [61]). 

We generate mock realisations in which UHECRs suffer simulated magnetic deflections 
d in the range 5° — 25° and determine the range of deflection angles for which a correlation 
with the source population can be established. We also compare these model distributions 
to the 69 observed PAO UHECRs. We simulate random magnetic deflections in the arrival 
directions of UHECRs by allowing for a shift in their positions with respect to the position 
of their sources. We do so using a 2 dimensional Gaussian function with a width d centred 
at the position of the source and according to it we generate randomly oriented angular dis- 
placements for each UHECR in each mock realisation. 

In figure 10 we show the distribution of values of X obtained for mock UHECRs, whose 
sources are nearby galaxies, and which have suffered magnetic deflections in the range 0°— 25°. 
We find that even if the sources of the UHECRs are correlated with nearby galaxies (as the 
ones in our model) increasing the deflection angle dilutes the correlation signal and for a 
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dataset of 69 events the correlation with the source population is wiped out by deflections 
greater than ~ 20°. We also show the value of X obtained for the 69 PAO observed UHE- 
CRs, XAuger- At present, given the value of XAuger obtained, we cannot exclude deflections 
of more than 20°. With a larger dataset, if the correlation of observed UHECRs with LSS 
persists, we will be able to exclude large deflections and hence a heavy composition. 

4 Discussion and Conclusions 

We have used the set of 69 observed PAO UHECRs with energy greater than 55 EeV to assess 
whether their arrival directions are correlated with nearby galaxies. We model the UHECR 
source distribution using the 6dF and PSCz galaxy surveys, taking into account the UHECR 
flux suppression due to energy losses during particle propagation and the non-uniform PAO 
exposure. We have shown that the probability that the observed UHECRs are consistent 
with an isotropic distribution is less than 5% for most of the parameter space considered. 

The analysis performed with the 6dF, which is being used here for the first time to 
model the local UHECR source distribution agrees well with that obtained with the PSCz 
despite the two surveys having different median depths and fields of view. The results ob- 
tained here are in agreement with those of other studies such as [24] in that the observed 
UHECRs exhibit an anisotropy which is consistent with the distribution of nearby galaxies, 
although our results are not directly comparable due to the different statistical approaches 
and UHECR source models used. 

In this work we have taken a Frequentist approach to the long standing question of 
the origin of UHECRs, as we consider it the most straightforward approach with the small- 
est number of assumptions regarding the source population. Other authors have taken a 
Bayesian approach to the question [30], and despite using different models of source popula- 
tion they have obtained results that are qualitatively similar to the ones obtained in this and 
similar analyses, namely that the UHECRs observed to date (or at least a significant fraction 
of them) are inconsistent with an isotropic distribution of sources and consistent with models 
in which the sources are luminous galaxies in the nearby universe. 

Our analysis has been performed with two galaxy surveys which have been carried out 
using different selection criteria and hence probe different astrophysical populations. Ob- 
serving a correlation with a specific astrophysical population does not constitute proof that 
the sources of UHECRs are members of that population, as matter in the universe is clus- 
tered and different astrophysical populations are correlated with each other. The confidence 
intervals with which we can exclude isotropy using the two surveys are comparable which 
exhibits the robustness of our results. With a larger UHECR dataset one may hope to dis- 
tinguish between different source populations, if the strength of EGMFs is not too large and 
the composition is nucleonic. 

We note that in this work we have only considered the case of faint sources which is 
well motivated in the absence of significant small scale clustering (repeaters) in the dataset. 

We have excluded the Galactic plane throughout our analysis, both due to the incom- 
pleteness of the galaxy surveys and because UHECRs travelling through the Galactic plane 
are likely to suffer significant deflections due to the Galactic disc's magnetic field. We have 
attempted to set a limit on the deflection suffered by the PAO UHECRs since the strength 
of Galactic magnetic fields and EGMFs as well as the composition of the UHECRs remain 
loosely constrained. We have shown that the arrival direction distribution of UHECRs which 
have been deflected by more than 20° would be difficult to distinguish from an isotropically 
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produced dataset with a dataset of ~ 70 events. With the present dataset we cannot exclude 
deflections above 20° but with a larger dataset it will be possible to do so if the observed 
correlation with LSS persists. Excluding large deflections would suggest the absence of Fe 
nuclei for most models of EGMFs. 
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